Phonon Effects on Spin-Charge Separation in One Dimension 
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Phonon effects on spin-charge separation in one dimension are investigated through the calculation 
of one-electron spectral functions in terms of the recently developed cluster perturbation theory 
together with an optimized phonon approach. It is found that the retardation effect due to the 
finiteness of phonon frequency suppresses the spin-charge separation and eventually makes it invisible 
in the spectral function. By comparing our results with experimental data of TTF-TCNQ, it is 
observed that the electron-phonon interaction must be taken into account when interpreting the 
ARPES data. 
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It is commonly accepted that most one-dimensional 
(ID) correlated electronic systems cannot be properly 
described by the traditional Fermi liquid theory, instead, 
their behaviors are well predicted by the Luttinger liq- 
uid theory 0, Q S 01- One of the key features of the 
Luttinger liquid is the spin-charge separation: the low- 
energy excitations are not quasiparticles with charge e 
and spin 1/2 together, rather, they are collective modes 
of spin and charge excitations separately, called spinous 
and holons. Since spinon and holon move with different 
speed, they eventually decouple. Following earlier works 
many studied spin-charge separation with var- 
ious theoretical schemes. In particular, one could explore 
the existence of the spin-charge separation by calculat- 
ing the spectral function 0, IE M S ^| j which has 
direct relation to the angle-resolved photoemission spec- 
troscopy (ARPES). In some recent works, the spinon and 
holon branches have been observed by ARPES performed 
on some ID materials such as S'rCuO2-0, 
Because both electron-electron and electron-phonon in- 
teractions exist in many low-dimensional materials, it is 
important to address the role of these interactions on the 
spin-charge separation. 0, El E E3, El 113 

The one-dimensional Holstein-Hubbard model (HHM) , 
which is the simplest model involving both electron- 
phonon (e-p) and electron-electron (e-e) interactions, has 
been used extensively to describe some low-dimensional 
materials. Since electrons in these materials are strongly 
correlated, the interplay between electron-phonon inter- 
action and Coulomb repulsion should have profound ef- 
fect on the spin-charge separation, and we expect to ob- 
serve these effects by investigating the single-particle ex- 
citation spectra. The spectral function provides valuable 
insights into the usually complicated many-body sys- 
tems, such as high-temperature superconductors, cuprate 
ladder compounds, and organic conductors. For example, 
very recently, by using the exact diagonalization method, 
Fehske et aZ.j2^ calculated the spectral function of the 
Holstein-Hubbard model on a finite system and found a 



Mott-insulator to Peierls-insulator transition at a com- 
patible ratio of the e-e to e-p interactions. 

Bearing these in mind, we compute the one-electron 
spectral function of the HHM by applying the rece ntly 
developed cluster perturbation theory (CPTllsL l2ll 221 
together with an optimized phonon approach. |2,'^ |2J, l25l | 
The spectral function at full frequency region with rich 
satellite structures is obtained in the model of both e-e 
and e-p interactions for the first time. Phonon effects 
on spin-charge separation are focused in the presence 
of e-e interactions from weak to strong coupling. It is 
found that the retardation effect due to the finiteness of 
phonon frequency does not favor the spin-charge separa- 
tion. In weak interaction regimes, a peak in the spectral 
function was observed which is consistent with the exis- 
tence of a metallic phase as proposed recently by Clay 
and Hardikarji^. Furthermore, it is observed that one 
must take electron-phonon interaction into account when 
interpreting the ARPES experimental data in the one di- 
mensional material. 

The HHM accounts for a tight-binding electron band, 
on-site Coulomb repulsion between electrons of opposite 
spin, and coupling of charge carriers to local phonons: 

H = y^(c| erCi+l^cr + H.C.) + t/ "-tT"a 

where c\ ^{ci^cr) creates (annihilates) an electron with a 
on site «, and b\ and bi are creation and annihilation oper- 
ators of the local phonon mode at site i, respectively, t is 
the electron hopping constant between nearest neighbor 
sites which will be set as the energy unit in our calcu- 
lations, LOf) is the bare phonon frequency, and g is the 
electron-phonon coupling constant. 

For the calculation of the spectral properties within 
the framework of CPT@,0|, one divides the lattice into 
fully equivalent clusters of a finite sites. For each cluster. 
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we calculate the Green's function Gij{z){= G^Az) 
G~j{z)), with G^j{z) defined as 
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where c. 



+ = J 



z±(H- Eo) 



10o), 



(2) 



Ci, and \(j>o) being the ground state 



of the cluster, which is obtained by using the Lanczos 
exact diagonalization (ED ) method within an optimized 
phonon approach [3 under open boundary condi- 
tions. Two terms in Gij corresponding to electron and 
hole propagation, respectively, can be obtained. The 
CPT treats the intercluster hopping by a strong-coupling 
perturbation, i.e., (t/U) expansion. The lowest-order 
CPT approximation to Green's function gives 



GcPT{k,z) = l^e-'=(-^)G,,,(iVfc,z), 

id 



(3) 



where Gij{Q, z) is the Green's function of the full system 
and N is the size of clusters. The spectral function is 
then A{k^uo) = —^Im[GcPT{k, z)], where z = uj + i-q 
with T] defines the width of peaks in the spectral function. 
Since the Fermi energy is set to zero, the spectral function 
has the symmetry of A{k,iS) = A{-k — k, —uj) due to the 
electron- hole symmetry of the model (1). 

In the absence of interactions, A{k, oj) obtained by the 
CPT method is exact |^ |^ , while for interacting mod- 
els, the accuracy of CPT depends on the size of the clus- 
ter and the number of optimal phonon chosen. To test 
the accuracy of the approach we use, we calculated the 
spectral function of the Hubbard model with exactly the 
same parameters as used by Benthien et al0 and ob- 
tained agreeable results. We also calculated the first two 
spectral momenta and they match exactly to those ob- 
tained by White |26|. Furthermore, based on our previous 
technique analysis j23, system parameters were carefully 
chosen in this work to ensure that our results mimic ther- 
modynamic limit. Results obtained in this Letter were 
for N = 6, r] = O.lt, and three optimized phonons at each 
site, with relative error 10~^ for the total energy. 

Three energy scales govern the physics of the HHM: 
the on-site Coulomb repulsion (U), the electron-phonon 
coupling (A = 2(7^/iu;o) and the bare phonon frequency 
(tJo)- The ground state of the system at half- filling is 
a Mott-Hubbard insulating (MI) state when U is large, 
and shows spin-density-wave (SDW) fluctuations. When 
electron-phonon interaction dominates, the system is in 
the Peierls insulating state (PI), characterized by the 
charge-density-wave where both spin gap and change gap 
are finite, while in the MI state, the spin gap vanishes. 
Very recently, it was reported that there is a metallic 
region intermediate between the PI and MI states with 
superconducting pairing correlation dominates. |2^ 

In the absence of phonons {g = 0), Eq. ^ is just 
the Hubbard model whose physics have been extensively 
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FIG. 1: The obtained spectral function A{k,u!) of (a) the 
large-ojo Holstein model and (b) the negative-U Hubbard 
model at half-filling. HO and SP stand for holon and spinon 
excitations respectively. 



studied and it is well known that in the Hubbard model, 
spin and charge separate. 0j S S On the other 
hand, when U = 0, Eq. ^ is another extensively studied 
model, the Holstein model (HM). In the strong electron- 
phonon coupling region, the ground state of the Holstein 
model at half filling is either a bipolaron insulating (BPI) 
state in the large Uq limit, characterized by the configu- 
ration where each site is either empty or doubly occupied 
because phonons produce an attraction between the elec- 
trons, or a traditional band insulating state in the small 
ujo limit. (20.. 28] In the weak coupling region, the Peierls 
gap is suppressed by the phonon quantum fluctuations 
and the ground state is at metallic (M) phase, [j^ E^ . Eflj | 
When the phonon degrees of freedom are integrated out, 
the spin-1/2 Holstein model could be mapped onto the 
Hubbard model with an effective dynamical attraction 
Ueff{ijj) = — A/(l — w^/cjq). Here one also expects to 
observe the spin-charge separation. In the anti-adiabatic 
limit, i.e., wp ~* oOi the attraction becomes instantaneous 
and equals to the bipolaron binding energy A, which has 
already been reported long time ago. ,2Sj Obviously, for 
any finite phonon frequency loq , one must consider the re- 
tardation effect fully which could not be simply presented 
by the above Ueff{uj). 

To have a sense of the magnitude of wq, we give a 
comparison of the spectral function of these two models 
in Figure 1. It is clearly shown that the spectral func- 
tion of the Holstein model is almost the same as that of 
the negative-U Hubbard model. Some minor differences 
due to finiteness of wq are invisible in the figure (in other 
words, tJo = 8 is almost at the antiadibatic limit). This 
result is not trivial as seen at first glance because it im- 
plies the single-particle excitation of the system with a 
large phonon frequency is similar to that at the antiadi- 
abatic limit, which is consistent with the existence of a 
quantum metal-insulator phase transition in the Holstein 
model. [Ulillsol The peaks labelled "SP" and "HO" in 
Fig. 1 refer to the spinon and the holon branches, respec- 
tively, signaturing the spin-charge separation. Compare 
to the conventional Luttinger liquids (e.g., the positive-U 
Hubbard model), the charge velocity (vp) is smaller than 
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FIG. 2: The spectral function ^4(0, u>) of the Holstein model. 



the spin velocity (wo-). 

Fig. 2 illustrates the retardation effects systematically. 
Starting from the strong coupling case (A = 2), we ob- 
serve that as we decrease wq from 4 to 2, the incoherent 
part of the spectra becomes more important. As a con- 
sequence, the spectral weights of the spinon and holon 
excitations are much smaller. This could be regarded 
as a continuation from Fig. 1: ujq — oo —>■ 8 —>■ A. 
When the phonon frequency is further reduced, the spec- 
tral weights correspond to phonon excitations become 
dominant, which is quite different from that in the an- 
tiadiabatic regime where the "spinon" and "holon" exci- 
tations are clearly the dominant ones. Therefore, due 
to the strong mixing of the coherent and incoherent 
excitations, it is difficult to single out the spinon and 
holon excitations, instead, one observes an almost flat 
band dispersion with exponentially small spectral weight. 
The dominant peaks in the incoherent part of the spec- 
tra are related to multiples of the (large) bare phonon 
frequency broadened by electronic excitations. Such 
electron-phonon mixed nature of excitations could be 
seen in the spectra away from the Fermi surface. 

It is quite natural to expect that the separation of spin 
and charge excitation will become smaller with the de- 
crease of the electron-phonon coupling strength. This is 
clearly reflected in the spectral function. As shown in 
the first row of Fig. 2, when we reduce the electron- 
phonon coupling A at Wo = 4, the difference between 
the spinon and the holon excitation at fc = becomes 
smaller and eventually invisible. As we reduce phonon 
frequency, retardation comes into play. It is also inter- 
esting to observe that there is an excitation split in the 
weak coupling case (Fig. 2, A = 0.5, ojq = 1). Such 
splitting is not due to spin-charge separation. In fact, 
by carefully comparing this spectral function with that 
of a spinless Holstein model at corresponding electron- 
phonon coupling, we found that the splitting is caused 
by the poiaron interaction. 

Figure 3 shows the spectra of up-spin electron for the 



(a) spin-1/2 il 

a/'\ 


(b) spinless A 
ffy 


nv 

ff( 1 






u\ 
















— — 








^r^^X^ 




— 

r-~^ 1 ~^ 1 . 


(0-|J -2 2 


(j3-|I -2 S 



FIG. 3: Spectral function A{k,ui) of the spinful (a) and spin- 
less (b) Holstein model. A = 0.5 and uio = 1. 
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FIG. 4: Spectral function ^1(0, ui) of the Holstein-Hubbard 
model at half filling, ujo ~ 1- 

spin-1/2 Holstein model at half- filling in comparison with 
the spinless Holstein model. In the weak coupling regime 
the two spectral are almost the same, indicating that the 
existence of the down-spin electrons have nearly zero ef- 
fect on the spectral function of the up-spin electrons. It 
shows that the phonon-mediated interaction between up- 
and down-spin electrons is very weak so they are almost 
decoupled in this case. Thus the splitting of the excita- 
tions cannot be attributed to the spin-charge separation, 
rather, it is due to the interaction among polarons. No- 
tice that there is a small peak labelled as "A" in the 
spectral weight, which is almost dispersionless in small k 
regime (see Fig. 3) and is suppressed by the on-site re- 
pulsion U (see Fig. 4). Since it is appeared in the metal- 
lic region intermediate between PI and MI phases, we 
speculate the peak "A" may be related with the electron 
pairing, as discussed recently by Clay and Hardikerp^. 
although further investigations are definitely deserved. 

Now we turn the Coulomb interaction U on and dis- 
cuss its effect on the spectral function in the pres- 
ence of electron-phonon interaction. From Figure 4, we 
also see that, the electron-electron interaction and the 
electron-phonon interaction have opposite effect on the 
spin-charge separation, as we intuitively expected. At 
given electron-electron interaction (for example, U = 
4), increasing the electron-phonon interaction tends to 
broaden the excitation bands and leave the spin-charge 
separation invisible. While at fix electron-phonon inter- 
action, the electron-electron interaction increases the sep- 
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FIG. 5: Density plot of the spectral function A{k,u!) for (a) 
the Holstein-Hubbard, and (b) the Hubbard model, (c) The 
spectral function A{0, uj) of corresponding models at different 
electron-phonon coupling. SH stands for the shadow band. 

LJo = 1. 



ration and found that the electron-electron interaction 
competes with the electron-phonon interaction on the 
spin-charge separation, and the retardation effect due to 
phonons may diminish the spin-charge separation in the 
spectral function. We also found polaron splitting and 
observed a peak that may related to electron pairing in 
the weak e-p coupling limit. 
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aration between the spin and the charge excitations. 

By further increasing the electron-electron interaction, 
we illustrate the role of electron-phonon coupling in 
Fig. 5. Comparing Fig. 5(b) with Fig. 5(a), one may 
still observe the signature of spin-charge separation as 
A increases, but the holon branch is much more broad- 
ened and its spectral weight decreases, while the spinon 
branch is nearly unchanged. This is clearly shown in 
Fig. 5(c). Consequently, the so-called shadow band, 
which orig inates from the diverging spin fluctuations at 
is nearly gone, due to the fact that the 
shadow band is actually the continuation of the holon 
band. 

Finally, we make a comparison on our results to the 
ARPES experiments of the quasi-one-dimensional or- 
ganic conductor TTF-TCNQ.[il[3 We notice that in 
the Fig. 7 of the Ref.^^, the charge branch (b) and 
the spin branch (a) are weakly connected with each 
other. The result can not be explained by the Hub- 
bard model alone(Fig. 5(b)). As our results suggest, the 
electron-phonon interaction must be taken into account. 
From Fig. 5(a), it is clearly seen that the dispersion of 
the charge branch is weakened in the spectral function 
and the spinon branch is weakly connected to the holon 
branch. A similar broadening due to phonon is also ob- 
served in one dimensional SrCu02.|2J| These facts indi- 
cate that one should expect significant contribution from 
the electron-phonon interaction to the spectra of these 
strongly correlated quasi-one-dimensional materials. Of 
course, our conclusions are based on numerical studies of 
Eq. (1) so to have detailed analyzes of experiments, one 
may use models different from Eq. (1) but we believe 
essential physics remain unchanged. 

In summary, by applying the CPT together with an 
optimized phonon approach, we have studied the spec- 
tral function of the one-dimensional Holstein-Hubbard 
model at half filling. A comprehensive picture for the 
spectral function in the presence of electron-electron in- 
teraction and electron-phonon interaction was presented. 
In particular, we addressed the issue of spin-charge sepa- 
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